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Abstract 

A formalism is developed which allows to determine the locations of all local symmetry axes 
of three-dimensional particles with overall icosahedral symmetry. It relies on the fact that the 
root system of the non-crystallographic Coxeter group ifs encodes the locations of the planes of 
reflection that generate the discrete rotational symmetries of the particles. Via an appropriate 
extension of the root system, new planes of reflection are introduced which determine local axes 
of rotational symmetry. An easy-to-implement formalism is derived that allows to compute the 
surface structure of any three-dimensional icosahedral particle with local symmetries. It can be 
used also for particles with overall octahedral and tetrahedral symmetry in conjunction with the 
root systems of the corresponding reflection groups. 

AppUcations to viruses are discussed explicitly. It is shown that the concept of quasi-equivalence 
in Caspar-Klug Theory corresponds to the special case of local six-fold symmetry axes contained 
in the theory developed here, and the corresponding geometries can hence be obtained with this 
formalism based on the root system of H^. 

Moreover, as a by-product, the theory answers the long-standing open question why only 
certain types of capsomeres, i.e. clusters of protein subunits, are observed in the surface structures 
of viruses. Since the types of the capsomeres are determined by the orders of the local symmetry 
axes on which they are located, the possible types of capsomeres are restricted by the spectrum 
of local symmetry axes allowed by the theory. Based on this we determine the spectrum of all 
capsomere types that may occur in viral capsids and give explicit examples for the lower-order 



1 Introduction 

Icosahedrally symmetric particles with local symmetries play a crucial role in virology, because viruses 
encapsulate their genome in protein containers, called viral capsids, that are organised in protein 
clusters centred on the symmetry axes of the capsids. These axes of symmetry may be global, i.e. 
corresponding to rotations that leave the whole structure invariant and are described by the overall 
group structure, or local, in which case only the local environments around the symmetry axes. An 
example of a local symmetry is the discrete rotation about the centre of the capsomere, that leaves 
only the capsomere itself, but not the whole capsid, invariant. The classification of all possible local 
symmetries of icosahedral surface structures hence determines the spectrum of possible viral capsids 
with icosahedral symmetry. 

Viruses with icosahedrally symmetric capsids occur predominantly in nature Q]. By making use 
of symmetry, a virus can efficiently encode the structure of its capsid: it only needs to encode the 
locations of the protein subunits modulo symmetry, that is only the locations of a subset of the 
proteins, and the locations of all other proteins follow by symmetry. Since the icosahedral group is 
the finite symmetry group of largest order in three dimensions (order 60, as opposed to order 48 for the 
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octahedral group, and order 24 for the tetrahedral group) it is hence not surprising that icosahedral 
symmetry plays a distinguished role.^ 

If only the symmetry axes of the icosahedral group (12 5-fold axes, 20 3-fold axes and 30 2- 
fold axes) could be used to encode the locations of protein subunits, only small viral capsids would 
occur. Viral capsids therefore also follow local symmetry axes. This has already been observed by 
Caspar and Klug 2 , and they classify the structure of particles that follow a particular type of local 
symmetry, which runs under the terminology quasi- equivalence in their work. However, experimental 
evidence has shown that viruses follow also local symmetries different from these fT^. Inspired by this, 
a mathematical formalism is developed here which can be used to determine the surface structures of 
all icosahedral particles with local symmetries. 

The approach adopted here is based on the fact that the rotational symmetries of the icosahedral 
group are generated by the reflections that form the non-crystallographic Coxeter group Hs. The 
locations of the corresponding planes of reflection are encoded by a set of vectors, called root vectors, 
that are the (normalized) vectors perpendicular to the reflection planes. In order to obtain further 
planes of reflection that correspond to local rotational symmetries and are compatible with the overall 
group structure (global symmetries), the set of root vectors has to be extended in an appropriate way. 
In order to achieve this, affine extensions of non-crystallographic Coxeter groups are used JHl as 
discussed in sectional They lead to an augmentation of the set of roots by a further root, and Z- 
linear combinations of this enlarged set of roots define nested point sets that are subsets of generalised 
lattices. The points in these sets correspond to the normal vectors of additional reflection planes 
and hence encode the locations of additional (local) symmetry axes that are compatible with the 
overall group structure. The structures of icosahedral particles with local symmetries can thus be 
reconstructed based on an affine extension of the non-crystallographic Coxeter group Hs. In section 
[21 a simple tool-kit is developed based on this formalism which can be used to compute the structure 
of these particles. It is easy to implement and can be used also by people who are not familiar with 
group theory. Applications are considered in sectional 

It is important to note that the formalism presented here can be used to construct all icosahedral 
particles with local symmetry axes. The reason for this lies in its connection with the projection 
formalism known from the study of quasicrystals JjQ, JJ and Penrose tilings In particular, 

it is known that point sets with local symmetries can be obtained via a projection from a higher 
dimensional regular lattice such as, for example, the root lattices corresponding to indecomposable 
finite root systems All local symmetries of particles with overall icosahedral symmetry can hence 
be obtained via projection from the root lattice of Dq JUj. Since the root system of Hs can be 
obtained from the root system of Dq by projection [T^, and the root lattice of Dq is given by Z-linear 
combinations of the root vectors, which after projection lead to Z-linear combination of the root 
vectors of Hs, one can equally work with Z-linear combinations of the root vectors of Hs themselves 
provided that one has a guide as to which linear combinations are important. This information is 
provided by the affine extension of the group, and leads to the tool-kit for the construction of particles 
with local symmetries presented in subsection 12.21 

Therefore, even though one could in principle construct the surface structures of the particles via 
projections from the root lattice of Dq^ a different approach is adopted here, which is less cumber- 
some because remnants of the higher dimensional space are only implicitly contained in the theory. 
All practical computations reduce to very simple operations in three dimensions that can easily be 
implemented on a computer. Applications are discussed in section |S1 

^It allows to compress the data necessary for encoding of the locations of the protein subunits in the capsid by a 
factor of 1/60. 
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2 The group theoretical set-up 



Due to the importance of icosahedral symmetry in virology, we present the theory with focus on this 
case. However, the theory can be used also for the derivation of the surface structures of particles 
with octahedral and tetrahedral symmetry, if the root systems of the corresponding reflection groups 
are used instead of the root system of H^. 

2.1 Background 

The Coxeter group Hs is the only non-crystallographic Coxeter group in three dimensions. It is a 
subgroup of i^4, the only non-crystallographic Coxeter group in four dimensions, and contains H2 as 
a subgroup |n|. The root systems of H4 and H2 can be obtained via projection from the root systems 
of the crystallographic groups A4 and Es^ respectively. As in the case of Hs and Dq^ the higher 
dimensional crystallographic root systems in these pairings are distinguished by the fact that all root 
vectors are of the same length. 

The elements of Hs are reflections rj at planes perpendicular to the root vectors (or roots in short) 
aj. The collection of all root vectors is called the root system, denoted as A, and encodes all elements 

in the group as follows: to each aj G A there corresponds a reflection rj{x) = x — ^ ^^^j^^j ^ a^, 
which is a reflection at the plane perpendicular to the root aj . The root system of Hs is given by the 
following 30 vectors 1^: 



equation x = x + 1. 

These vectors point from the centre of an icosidodecahedron to its vertices, as shown in Fig. Q on 
the left, and encode each the location of a plane of reflection, as shown on the right. The intersections 



Figure 1: Examples of planes of reflection encoded by the root vectors of Hs (left) and the root 
polytope that encodes the locations of the planes of reflection (right). 

of the planes of reflection mark the locations of the axes of rotational symmetry of the icosahedral 
group. For example, the intersection of the two planes in the figure corresponds to two axes of 5-fold 
symmetry, which intersect the sphere at two of the 12 5-fold vertices of the (spherical) icosahedron. 
The other reflection planes are not shown, but their intersections with the surface of the sphere are 
indicated as geodesies (spherical arcs obtained as intersections of planes through the origin with the 
surface of the sphere). The intersections of the geodesies mark the locations of all symmetry axes, 
and one can hence reconstruct the locations of all 12 five- fold, 30 two- fold and 20 three- fold axes from 
them. 
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In order to extend the root system in a way compatible with overah icosahedral symmetry and 
hence find new planes of reflection as well as the associated axes of local symmetry, one needs to 
observe that the information encoded by the root system A can further be condensed into three so- 
called simple roots. They form a basis from which all other root vectors can be constructed by positive 
or negative Z[r]-linear combinations, where Z[r] = {a-\-rb\a^b G Z} corresponds to an extended ring 
of integers. 

In particular, a possible choice of simple roots in the orthonormal basis is 

ai = (0, 0, 1) , a2 = |(-r^ -r, -1) , as = (0, 1, 0) . (2) 



The relations between the simple roots are encoded in the Cartan matrix C, 

V Kki) \o _r 2 I 



Via an afiine extension this matrix is extended by an additional row and column, which encode 
a further root, that corresponds to an affine reflection. It has been shown in ^Hl? that the affine 
extended Cartan matrix is given by 
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and that the affine reflection act as a translation T by the highest root an = Tai + 2ra2 
—t'(j02 = (1,0,0). The three other reflections are cyclic operations of order two and the products of 
any two of them correspond to rotations around the origin: 
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(5) 



The affine extended group is hence generated by the three reflections ri, r2, rs as well as T. 

Three-dimensional points sets with local symmetries can be constructed via an iterated action of 
the generators of the extended group on the origin. If the action of the translation operator T is not 
restricted, space is densely filled in this way. However, if T acts only a finite number of times, N say, 
while the action of all other operations is not restricted, point sets S{N) are obtained which are finite 
subsets of cut-and-project quasicrystals with centrally symmetric acceptance windows. For increasing 

the patches become larger and more dense, and correspond to cut-and-project quasicrystals with 
increasingly larger acceptance windows. Examples have been worked out explicitly for the two- 
dimensional subgroup H2 in 18 , and for in ^^l^ where applications to carbon cage structures 
have been considered. 

Due to the fact that that T acts as a translation by the highest root, a simple expression for these 
point sets can be given in terms of the root system A (see for the case of Hs): 



'= \ ^n,a|n, gNU{0}, 



< N 



(6) 



i.e. S{N) consists of all points that are obtained as an N U {0}-linear combination of up to N roots 
in A. These sets form the starting point for the construction of icosahedral particles with local 
symmetries in the next subsection. 
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2.2 The tool-kit 



The sets S{N) in decompose into disjoint subsets Sk{N) of vectors of equal length that are 
invariant under the action of Hs: 

SiN) = U^-,Sk{N), (7) 

where Kn denotes the number of different subsets of vectors of equal length in S{N). These subsets 
are of different cardinalities, and a given set Sk{N) may increase in cardinality for increasing 
i.e. one has the set inclusion Sk{N) ^ Sk{N -\- 1) for / G N. 

Each set Sk {N) encodes the surface structure of a particle that has local symmetry axes compatible 
with the overall group structure. Since only the orientations but not the length of the vectors in the 
set is important for the reconstruction of the locations of the local reflection planes, two sets Sk{N) 
correspond to the same particle if they contain the same vectors modulo length. Therefore the range 
of different (modulo length) sets Sk{N) corresponds to the spectrum of different icosahedral particles 
with local symmetries. The tool-kit given here provides a method for the systematic evaluation of 
this information and hence for the construction of the surface structures of all such particles. 

By definition, 5(1) corresponds to the root system A. The reflection planes orthogonal to the roots 
intersect precisely at the symmetry axes corresponding to the (global) symmetries of the icosahedral 
group, and it is hence a structure without local symmetries. Therefore, local symmetries are encoded 
in the point sets S{N) with N > 2. For each the sets Sk{N) in {TJ need to be evaluated individually. 
This is done by an inductive process starting with the case N = 2. It is presented in the form of 
a tool-kit which facilitates easy computational implementation and can be used also by people not 
familiar with the group theoretical set-up: 

1. Add any two vectors in Q and group these vectors into sets <S/c(2), /c = 1, . . . , i^2, of vectors of 
equal length. For each /c, compute the planes through the origin that are perpendicular to the 
vectors in Sk{2) according to the equation 

3 

^VjXj=0 (8) 

for any vector v G Sk{2) with Cartesian components vj. The intersections of these planes 
determine the locations and types of the corresponding local symmetry axes. They may either 
be computed analytically or graphically displayed with a visualisation tool such as Maple. 

2. Compute the structure of S{N) inductively based on the structure of Sk{N — 1) as follows: Add 
any N vectors in (^) and group them by their length into sets (N). If the length of the vectors 
in a set {N) corresponds to that of one of a sets at iteration step N — 1^ Sk{N — 1) say, form a 
new set that contains both sets, i.e. introduce a new set Sk{N) := {N)[JSk{N — 1). Otherwise, 
add the set as a new set to the list, i.e. Si{N) = (N)^ where the index / is chosen such that 
sets are numbered consecutively. The collection of sets obtained in this way corresponds to the 
sets Sk{N)^ k = l,...,i^Ar in 0^. For each of them, compute the corresponding symmetry 
planes according to as before. 

Applications of this tool-kit are discussed in section |31 Note that the intersections of the reflection 
planes with the surface of the sphere lead to spherical tessellations of the surface of the sphere. The 
set inclusions of the type Sk{Ni) C Sk{N2) for N2 > Ni imply that the spherical tiling (tessellations) 
corresponding to these sets are related by self similarity: larger copies of the tiles (areas) implied by 
the spherical tiling related to Sk{Ni) are replaced by an integer number of rescaled versions thereof in 

^As mentioned before, two sets of vectors that are identical modulo length represent identical particles. However, 
none of them can be discarded from the list as they may lead, at a later iteration step, to different point sets when 
further points are concatenated to them. 
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the spherical tihng related to Sk{N2). An example is shown in Fig. [21 where each of the four triangles 
in the original triangulation (with vertices shown as blue dots) has been replaced by four rescaled 
(smaller) versions in the refined triangulation. 

The tiling is shown only in one of the 20 faces of the spherical icosahedron (or, equivalently, in 
three copies of the fundamental domain of the icosahedral group), because the remainder of the tiling 
follows from this via the actions of the icosahedral group. Planes represented by lines connecting the 
blue dots arise from the vectors in the set Sk {Ni ) , while planes represented by dashed lines connecting 
the red dots correspond to the vectors that are contained in Sk{N2) in addition. 




Figure 2: Three axes of local symmetry shown in blue, additional symmetry axes added to the set for 
larger N are shown in red. 

This phenomenon is also known from the study of quasicrystals ^5]. Since quasicrystals are ob- 
tained from a projection formalism that is related to this construction as discussed in the introduction, 
this is not surprising. 



3 Applications to viruses 

In this section it is demonstrated how the theory can be applied in virology. In subsection 13.11 the 
tool-kit is applied to derive the surface structures of new viral particles with local symmetries. In 
subsection 13.21 it is shown that the theory answers, as a by-product, an open question in virology: it 
explains why only certain types of capsomeres (protein clusters) are observed in viral capsids. 



3.1 New particles with local symmetries 

The local symmetries of particles with icosahedral symmetry can be obtained from the root system 
of the non-crystallographic Coxeter group Hs via the tool-kit presented in the previous section. Ac- 
cording to this, we start by considering the set S{2) of all vectors obtained via a linear combinations 
of any two root vectors in and decompose this set into subsets <S/c(2) of roots of equal length. The 
types of the local symmetry axes that are obtained as intersections of the reflection planes encoded 
by the vectors in a set <S/c(2) depend on the relative orientations of these vectors, as can be seen from 
the examples discussed below. 



1. A particle with local 3-fold axes: 

Consider the subset <S/c(2) of S{2) given by linear combinations of any two vectors pointing to 
vertices of the icosidodecahedron with a joint edge, such as, for example, the vertices shown in 
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blue in Fig. [31 The resulting vectors intersect the edges at the points shown in red. In order 
to reconstruct the locations of the reflection planes orthogonal to these vectors, note that the 
organisation of the edges reflects the subgroup structure of they are grouped in decagonal 
arrangements according to copies of the subgroup H2. Each such decagon encodes 10 root 
vectors of , which in turn encode the symmetry planes that intersect at a global 5-fold vertex 
and its antipode on the sphere. Therefore, the bisection of the edges in the root vector picture 
in Fig. Olon the left correspond to symmetry planes that bisect the angles between two planes 
meeting at a 5-fold vertex at an angle 27r/10. The geodesies corresponding to the intersection 
of such planes with the sphere are shown in Fig. |31on the right superimposed on the spherical 
icosahedron. In order to keep the presentation transparent, only the geodesies in one of the 
20 faces of the spherical icosahedron are shown: geodesies arising from bisections of the angles 
within this spherical triangle are shown as dashed lines; moreover, the geodesies arising from 
bisections of angles in different spherical triangles (faces) are shown as dotted lines in order to 
show their impact on the spherical triangle under consideration. 



Figure 3: Linear combinations in the root system (left) leading to a spherical tiling with local 6-fold 
symmetry axes at the green vertices (right). 

Note that the intersections of these geodesies determine the locations of the local symmetry 
axes: there are three local symmetry axes of order six in the spherical triangle. Since this 
spherical triangle contains three copies of the fundamental domain of the icosahedral group, this 
setting corresponds to one new local 6-fold symmetry axis that is reproduced 60 times under the 
action of the icosahedral group. Note that the planes of symmetry, respectively the geodesies 
corresponding to them, not only fix the locations of the symmetry axes but also the relative 
orientation of the planes of symmetry and hence the surface structure of the corresponding viral 
capsids is completely determined by this formalism. 

2. A particle with local 5- fold axes: 

In order to obtain a particle with local 5-fold symmetry axes, a different subset Sk(2) of S{2) 
has to be considered, which inherits its order of symmetry from the global 5-fold axes. For this, 
consider linear combinations of root vectors that are located on two different i^2-orbits, as for 
example the vectors pointing to the vertices shown in blue in Fig. ^ on the left. These vertices 
are located on two i^2-orbits that are shown in green and red, respectively, and intersect at the 
vertex shown in yellow. Fig. ^ on the right shows how this picture translates into reflection 
planes. In particular, part of the geodesies corresponding to the reflection planes perpendicular 
to the 10 root vectors pointing to vertices on the green, respectively red, i^2-orbit are shown 
schematically in green, respectively red, around 5-fold vertices in the spherical tiling. The 
geodesies corresponding to the two root vectors pointing to vertices marked by blue dots in the 
flgure on the left are shown in blue in the spherical tiling on the right. Fig. Whence provides a 
dictionary on how to translate the geometry of the root vectors into geodesies on the sphere. 




7 



Figure 4: A dictionary for translating the i^2-orbit structure in the root system (left) into geodesies 
corresponding to reflection planes (right). 

Based on this, it is possible to construct the axes of local five- fold symmetry as follows. The linear 
combinations of any two root vectors corresponding to the blue dots in Fig. [SI (left) corresponds 
to a vector that points from the centre of the icosidodecahedron to any of the vertices marked 
in red in Fig. (left). Reconstructing the corresponding planes of symmetry as before leads to 
a network of geodesies with intersections of order 10 rather than 5. However, in order to obtain 
the spherical tiling that represents the surface structure of the viral particle, note that also the 
global 5-fold axes appear as 10-fold axes in this representation, because the angles between the 
five refiection planes at each 5-fold vertex are bisected by the refiection planes originating at 
the 5-fold vertex at its antipode on the sphere. Similarly, the local symmetry axes have to be 
interpreted as 5-fold symmetry axes in disguise, and every second geodesic has to be ignored. 
This leads to two different spherical tilings, shown in Fig. in the middle and on the right, 
respectively. 




Figure 5: Linear combinations in the root system (left) leading to two spherical tilings with local 
5-fold symmetry axes at the green vertices (middle and right). 

In order to visualise the structure of the viral capsids encoded by these tilings, the rules of 
the tiling approach need to be used [221 • Protein subunits are located in the corners of the 
tiles that meet at the local symmetry axes and are represented schematically as dots. The 
colours of the dots encode which proteins are mapped onto each other under the action of the 
overall symmetry group; the proteins themselves are identical and hence indistinguishable. Tiles 
with two and three dots represent, respectively, dimer- and trimer-interactions, i.e. interactions 
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between, respectively, two and three protein subunits. They are visualised schematically as 
spiral arms and may correspond for example to C- terminal arm extensions of the proteins. 
With this interpretation, the tiling in Fig. [HI on the right leads to the surface structure shown 
in Fig. El 



Figure 6: Schematic representation of the surface structures encoded by the tiling in Fig. (right). 

Note that the predictions of this tiling differ both by the orientations of the protein clusters 
(pentamers), which are rotated by an angle of 27r/5 with respect to the other alternative, and 
by the types and locations of the intersubunit bonds that stabilize the capsid. The tiling in the 
middle of Fig. O corresponds to the tiling in 20 and represents viral particles in the family of 
Papovaviridae. Viral particles following the alternative in Fig. are not known at present. 

Note that these sets are only two examples of sets in the decomposition of S{2). They have been 
chosen because they are the most instructive examples. There exists, for example, also a subset that 
does not introduce any local symmetry axes: it corresponds to linear combinations of each root vector 
with itself. These vectors point to the vertices of a larger copy of the icosidodecahedron, and the set 
hence encodes only the global symmetry axes (like the root system itself) and hence contains no new 
information. 

While the geometry of the particle in item 1 is known already from the Caspar-Klug classification 
(it belongs to a T = capsid), the geometries in item 2 are new. They correspond to particles 
with only 5-fold axes of symmetry (global and local). Due to overall icosahedral symmetry, any local 
5-fold axis occurs in multiples of 60, so that they are indeed the smallest geometries corresponding to 
particles with local five- fold axes. A classification of all icosahedral particles with local symmetries 
based on the theory developed here is in progress [2^ . 

3.2 Orders of the capsomeres 

The centres of the capsomeres in viral capsids correspond to the local and global axes of symmetry. 
Therefore, the spectrum of different capsomeres that may occur in nature is determined completely 
by the orders of the symmetry axes that can occur from a mathematical point of view. The group 
theoretical formalism developed here implies that the orders of all local symmetry axes must corre- 
spond to multiples of those of the global ones, because they inherit their symmetry from the structure 
of the root system due to the fact that they arise from Z-linear combinations of the root vectors. 
Therefore, the orders of the local symmetry axes must be Z- multiples of 2, 3 or 5, and local symmetry 
axes of, for example, order 7 cannot occur. Moreover, since smaller particles are obtained from linear 
combinations of only a few root vectors, multiples of more than twice the order of the global axes do 
not occur for them. Hence, for all smaller icosahedrally symmetric particles, which form the majority 
of viruses observed to date, one would expect a spectrum of local symmetries of order 3, 5, 6 and 10. 
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All of these are known to occur: order three has been observed for bacteriophage MS2 j^Sl, order five 
for polyomavirus ^S] and Simian Virus 40 JT], order six in all Caspar-Klug cases, and order 10 for 
L-A virus of yeast Pj. All these cases follow the construction principle outlined above. 

4 Discussion 

The formalism presented here is used in [22 to describe the surface structures of viral capsids in terms 
of tilings. The vertex sets of the tilings contain the points, that correspond to the local and local 
symmetry axes, and mark the locations of the capsomeres. The edges of the tiles follow the geodesies 
that arise as intersections of the planes of (local and global) symmetry with the surface of the sphere. 
In this sense, the tiling picture and the group theoretical picture in terms of local symmetry axes are 
equivalent ways of modelling the surface structures of viral capsids. 

This may suggest that the tiles are of secondary importance and that the symmetry axes should 
be the focus of attention. However, in this way, important additional information on the structure 
of viral capsids would be neglected as tiles represent the interactions between subunits belonging 
to different capsomeres and are hence biologically meaningful objects. In particular, the measures 
and shapes of the tiles are correlated with the strengths and types of the intersubunit bonds they 
represent [23 . Their structure is hence very important information in order to determine the types 
of the interactions (i.e. whether they are trimer- or dimer-interactions, that is interactions between 
three, respectively two, protein subunits), as well as their exact locations. This information is crucial 
input for the construction of assembly models [3 IHj i that investigate how viral capsids are assembled 
from their building blocks, as well as for the modelling and classification of additional structures such 
as the covalent bonds responsible for crosslinking [23| . 
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